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Next-to-leading order calculations of heavy quark three-jet production in e + e — 
annihilation are reviewed. Their application for the measurement of the b-quark 
mass at LEP /SLC and to the test of the flavour independence of the strong coupling 
constant are discussed. Prospects for future improvements are studied. 



1 Introduction 

Effects of the bottom-quark mass, m*, were already noticed in the early testsB 
of the flavour independence of the strong coupling constant, a s , in e + e~- 
annihilation at the Z-peak. Motivated by the remarkable sensitivity of the 
three-jet observables to the value of the quark mass, the possibility of the 
determination of mi, at LEP, assuming universality of the strong interactions, 
was consideredp. This question was analyzed in detail inp, where the necessity 
of the next- to- leading order (NLO) calculation for the measurements of mi, 
was also emphasized. 

The NLO calculation for the process e + e~ — * 3jets, with comnlete quark 
mass effects, has been performed independently by three groups ffiirO. These 
predictions are in agreement with each other and were successfully used in the 
measurements of the 6-quark mass far above threshpld&lj and in the precision 
tests of the universality of the strong interaction Bum at the Z-pole. In this 
talk we make a short review of such calculations. Furthermore, prospects for 
future improvements are discussed. 

It is surprising that at high energies the bottom-quark mass could be 
relevant, since it appears screened by the center of mass energy, m^/m^ — 10~ 3 
at LEP. Nevertheless, when more exclusive processes than a total cross section 
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are considered, like a n-jet cross section, mass effects can be enhanced as 
m 2 /m z /y c , where y c is the parameter that defines the jet multiplicity. 

Since quarks are not free particles and, therefore, their mass can be con- 
sidered like a coupling constant, one has the freedom to use different quark 
mass definitions, e.g. the perturbative pole mass Mf, or the MS scheme run- 
ning mass mb(n). Physics should be independent of it but at a fixed order in 
perturbation theory there is a significant dependence on which mass definition 
is used, as well as on the renormalization scale fi. The inclusion of higher 
orders to reduce these two uncertainties, due to mass definition and /i scale, is 
mandatory for an accurate description of mass effects. 



2 Three jet observables and the measurement of m& 

The observable proposed somej-time ago to measure the bottom-quark mass at 
the Z- resonance was the ratio El 



^3 = w — ,w > I 1 ) 



where T^- and T q are the three-jet and the total decay widths of the Z-boson 
into a quark pair of flavour q in a given jet-clustering algorithm. More precisely, 
the measured quantity is 

pM _ r 3j(z/c)/r b q s (/i) / a s {ij) \ 

^3 = TT1 VFFi = 1+ ao(y c )+rfe &o(r6,y c ) + bi{r b ,y c ) \ , (2) 

r3j(j/c)/r £ 7T v it j 

where now the sum of the contributions of the three light flavours i — u,d, s 
is included in the denominator. The R^ d and the R^ e observables differ only 
by the function ao (which is zero for R^ d ). This contribution originates from 
the triangle diagrams!^. It is numerically very small (0.002 for the Durham 
jet algorithm) and almost independent of the 6-quark mass. The bo and b\ 
functions give respectively the leading-order (LO) and NLO mass corrections, 
once the leading dependence on = M^/m 2 z , where M\, is the bottom-quark 
pole mass, has been factorized out. 

Ratios of differential two-jet rates, where the two-jet width T23 is calculated 
from the three- and the foutifit fractions through the identity T23 = V — T^j — 
T^j, have been studied intL3'Q. Ratios of event shape distributions have also 
been considered 0. 

Using the known relationship between the pole mass and the MS scheme 
running mass, 



M£ = mt{n) 



(3) 
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Figure 1: The observable Rg as a function of y c at the NLO for the four algorithms consid- 
ered. The dotted lines give the observable computed at the NLO in terms of the pole mass 
Mj, = 4.6 GeV while the dashed lines are obtained when it is written in terms of the running 
mass m^irnz) = 2.83 GeV. In both cases the renormalization scale is fixed at fi = niz, and 
cts{mz) = 0.118. For comparison we also plot in solid lines the LO results for Mj> = 4.6 GeV 
(LO-Mf,) and m b (m z ) = 2.83 GeV (LO-m 6 (m z )). 



we can re-express Eq. (Q) in terms of the running mass mfe(/x). Then, keeping 
only terms of order 0{a s ) we obtain 

< = l + ^aQ(y c ) + f 6 (/x)f&o(r6,yc) + ^6i(f6,yc,M)) , (4) 

7T \ TT / 

where fb{fi) = m^(/i)/m| and b\ = b\ + 2fo (4/3 - logr b + log(^ 2 /m|V). Al- 
though at the perturbative level both expressions, Eq. (@) and Eq. ([!]), are 
equivalent they give different answers since different higher order contribu- 
tions have been neglected. The spread of the results gives an estimate of the 
size of higher order corrections. 
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In fig. [l] we present our results for R\ in the four clustering algorithms: 
EM ^ Jade, E and Durham. For all the algorithms we plot the NLO results 
written either in terms oftll the pole mass, M& = 4.6 GeV, or in terms of the 
running mass at mz, mbimz) = 2.83 GeV. The renormalization scale is fixed 
to fj, = mz and a s (mz) — 0.118. For comparison we also show i?g at LO 
when the value of the pole mass, Mb, or the running mass at mz, mb(rnz), is 
used for the quark mass. 

Note the different behaviour of the different algorithms. In particular the 
E algorithm. As already discussed in □, in this algorithm the shift in the 
resolution parameter produced by the quark mass makes the mass corrections 
positive while by kinematical arguments one would expect a negative effect, 
since massive quarks radiate less gluons than massless quarks. Furthermore, 
the NLO corrections are very large in the E algorithm and strongly dependent 
on y c . All this probably indicates that it is difficult to give an accurate QCD 
prediction for it. For the Jade algorithm the NLO correction written in terms 
of the pole mass starts to be large for y c < 0.02. Note, however that the NLO 
correction written in terms of the running mass is still kept in a reasonable 
range in this region. Durham, in contrast, is the algorithm that presents a 
better behaviour for relatively low values of y c while keeping NLO corrections 
in a reasonable range. 

The theoretical predictions for the observables studied contain a residual 
dependence on the renormalization scale fi. To give an idea of the uncertainties 
introduced by this we plot in fig. ^a the observable R\ l as a function of /i for a 
fixed value of y c . Here we only present plots for the Durham algorithm, the one 
with the better behaviour. We use the following one-loop evolution equations 

a(M) = , m b (jM) = m h (mz) K~ l0 ^° , (5) 

A 

where a(/x) = a s (fi)/n, K = 1 + a(m^)/3 log(^ 2 /m|) with /3 = (11 — 
2/3Nf)/4, 70 = 1 and Np = 5 the number of active flavours, to connect 
the running parameters at different scales. 

Conversely, for a given value of R^ e we can solve Eq. (||) (or Eq. (Q)) with 
respect to the quark mass. The result, shown in fig. |b for Rf{y c = 0.02) = 
0.973, depends on which equation was used and has a residual dependence on 
the renormalization scale fj,. The curves in fig. are obtained in the following 
way: first from Eq. (|4|) we directly obtain for an arbitrary value of /x between 
mz and m^/10 a value for the bottom-quark running mass at that scale, 
mh(/Lt), and then using Eq. (|5|) we get a value for it at the Z-scale, mb(mz)- 
Second, using Eq. (|^) we extract, also for an arbitrary value of \i between mz 

b A modification of the standard Jade scheme, convenient for massive parton calculations!^. 
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Figure 2: Durham algorithm: (a) Renormalization scale dependence of the NLO predictions 
written in terms of either the running mass, NLO-mjfm^), or the pole mass, NLO-Af;,, for 
the ratio at a fixed value of y c . (b) Predicted value of mb(mz) for a fixed value of 
using either the running expression, NLO-mjfm^), or the pole expression, NLO-Af;,, as a 
function of the scale /i. 



and mz/10, a value for the pole mass, Mb- Then we use Eq. (||) at fi = Mb and 
again Eq. (^|) to perform the evolution from fi — Mb to /z = raz and finally 
get a value for mb(mz)- The two procedures give a different answer since 
different higher orders have been neglected in the intermediate steeps. The 
maximum spread of the two results, in this case of the order of ±200 MeV, 
can be interpreted as an estimate of the size of higher order corrections, that 
is, of the theoretical error in the extraction of mb(mz) from the experimental 
measurement of R$. 

These calculations were used by the DELPHI Coll. i to extract m,b{mz) 
from the experimental measurement of R$ e , see fig. ||[ and the result was inter- 
preted as the first experimental evidence (at 2-3 sigmas) for the running of a 
fermion mass since the data, are better described by the NLO-m^mz) curve. 
Also recently the SLD Coll. El has presented results for mb{mz)- The SLD anal- 
ysis is compatible with the DELPHI measurement. Nevertheless, the central 
values of m,b(mz) obtained from different clustering algorithms are scattered in 
the range Amb(mz) = ±0.49 GeV. This is probably due to the fact that E-likc 
algorithms, that are mainly used in this analysis, have huge NLO corrections 
thus making accurate QCD predictions difficult. 



5 




Figure 3: Experimental measurement of R^ e by the DELPHI Coll. in the Durham algorithm 
and preliminary results in the Cambridge algorithm compared with our NLO calculation 
written in terms of the running mass at the Z peak, NLO-mjfm^), or in terms of the pole 
mass, NLO-M5. The LO predictions are also plotted. 



3 Tests of the flavour independence of the strong interaction 

Assuming a given 6-quark mass the NLO calculation of heavy quark three-jet 
production cross section can be used to perform an improved test of the flavour 
independence of the strong coupling-constant. Such analysis was done for the 
first time by the DELPHI Coll. inn by using the R^ e observable defined in 
the Durham algorithm. Recently, the OPAL Coll. has presented a similar 
analysis by using different ratios of event shapes,distributions: D2, 1 — T, Mh, 
Bt, Bw and C. Instead, SLD has presented El results by analyzing the R3 
ratio in the E, EO, P, PO, Durham and Geneva algorithms. All these results are 
consistent with unity. No flavour dependence has been observed. Furthermore, 
the inclusion of mass effects is mandatory to achieve such agreement. 

4 Improving the &-quark mass measurements: the Cambridge al- 
gorithm 

The Cambridge algorithm has been introduced very recently^ in order to reduce 
the formation of spurious jets formed with low transverse momentum parti- 
cles that appear in the Durham algorithm at low y c . Therefore, compared to 
Durham, it allows to test smaller values of y c while still keeping NLO correc- 
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tions relatively small. Both algorithms are defined by the same recombination 
procedure and the same test variable 

y l3 = 2min(Ef, E])(l - cos%)/s , (6) 

where Ei and Ej denote the energies of particles i and j and 0y is the angle 
between their momenta. The new ingredient of the Cambridge algorithm is 
the so called ordering variable 

Vij = 2(1 - cos%) . (7) 

In the Cambridge algorithm one first finds the minimal t>y and then tests yij . 
If Vij < Vc-, the i and j particles are recombined into a new pseudoparticlc of 
momentum pk = pi + Pj but if > y c , the softer particle is resolved as a jet. 
The net effect of the new definition is that NLO corrections to the three-jet 
fraction become smaller. 

In fig. | we present the preliminary results from the DELPHI Coll.Ej for the 
R^ £ raticuiefined in the Cambridge algorithm and compare it with our NLO cal- 
culationEj written in terms of the running mass at the Z peak, NLO-mb(mz), 
or in terms of the pole mass, NLO-Mj,, for fi — vnz and a s (mz) = 0.118. As 
in Durham, the NLO-mb(mz) gives the best agreement. Furthermore, data 
are still compatible with the LO-m^mz) showing that the bulk of higher or- 
der corrections is described by the running of the b-quark mass. In contrast, 
although data could also be well described by the NLO-Mb curve, the NLO 
corrections become large when the pole mass paxametrization is used. 

The studies of the NLO-rab(mz) curve showlia that it is remarkably stable 
with respect to the variation of the scale fi. For the range m^/10 < fi < mz 
the estimate of the error in the extracted mb{mz) is reduced to ±50 MeV in the 
Cambridge scheme, with respect to ±200 MeV for Durham (±125 MeV if only 
NLO-mf,(m^) is considered). In contrast, when the NLO-Af;, parametrization 
is used we get ±240 MeV but strongly dependent on the lower n used. 

5 Conclusions 

In the last few years an important progress was done in the description of 
the Z-boson decay into three-jets with massive quarks. Next-to-leading order 
calculations have been done by three groups and have been successfully used in 
the analysis of the LEP and SLC data where mass effects have been clearly seen. 
Further studies of different observables and different jet-algorithms are oriented 
on the reduction of the theoretical uncertainty. One good candidate might 
be the Cambridge jet-algorithm, where the NLO corrections are particularly 
small and where the predictions in terms of the running mass, mi,(mz) are 
particularly stable with respect to the variation of the renormalization scale. 
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